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Temperature dependence of the Ua{ 1) anomaly is investigated by taking into account mesonic 
fluctuations in the U( 3) x U{ 3) linear sigma model. A field dependent anomaly coefficient function 
of the effective potential is calculated within the finite temperature functional renormalization group 
approach. The applied approximation scheme is a generalization of the chiral invariant expansion 
technique developed in [G. Fejos, Phys. Rev. D 90 , 096011 (2014)]. We provide an analytic 
expression and also numerical evidence that depending on the relationship between the two quartic 
couplings, mesonic fluctuations can either strengthen or weaken the anomaly as a function of the 
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of the chiral expansion is also discussed in detail. 

PACS numbers: 11.30.Qc, 11.30.Rd 

Keywords: Axial anomaly, chiral symmetry breaking, functional renormalization group 


I. INTRODUCTION 

Axial anomaly is the anomalous breaking of the Ua{ 1) 
subgroup of approximate f/z,(/V/) x Un{Nf) chiral sym¬ 
metry of quantum chromodynamics. It can be under¬ 
stood theoretically through instanton solutions of the 
classical equations of motion describing vacuum to vac¬ 
uum amplitudes with different topological winding num¬ 
bers [ji], 0. Although the origin of the anomaly has been 
clarified for a long time, very little is known about its 
finite temperature restoration, especially around and be¬ 
low the critical point. It is well established that due 
to disappearing instanton density at high enough tem¬ 
perature or chemical potential the 1 / 4 ( 1 ) anomaly has to 
vanish 0 0], however, it is an open question how relevant 
it is at the chiral transition point. There are recent ex¬ 
perimental findings that show a reduction of the rf mass 
near the chiral crossover temperature 0, which might be 
related to the restoration of the Ua{ 1) factor around Tc- 

The relevance of the anomaly around the critical point 
is crucial from the point of view of the order of the chiral 
transition itself, since the effective Nf flavor low energy 
description provided by the U(Nf) x U(Nf) linear sigma 
model 00 predicts a first order transition for Nf > 2, 
if the anomaly is restored at Tc and no explicit symme¬ 
try breaking terms are present. This is based on a fixed 
point analysis of the renormalization group (RG) flows 
Q, and also on explicit calculations of the effective po¬ 
tential 00. Even though there are indications that 
the RG argument might not survive for Nf = 2 EH3, 
it seems to remain correct for Nf > 3, which includes the 
most important Nf = 3 case. 

Beyond the anomaly free fixed point analysis of the 
aforementioned U(Nf) x U(Nf) scalar model, it can be 
easily extended with a 1 / 4 ( 1 ) breaking term, also known 
as the’t Hooft determinant. Its coupling is usually cho- 
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sen to be proportional to the topological susceptibility 
(Xtop ), in order to reproduce the Witten-Veneziano re¬ 
lation (WVR) [bil . Il5| . and hence shows an exponential 
suppression as a function of the temperature. It is based 
on linking Xtop to the instanton density at finite tem¬ 
perature, which, if the tunnelling amplitude is estimated 
by a semiclassical approximation, shows an exponential 
damping. However, this is controversial since strictly 
speaking, the semiclassical approximation is applicable 
only for temperatures higher than Tc- There have been 
attempts to use a modfied WVR 0 , but if the temper¬ 
ature dependence of the ‘t Hooft coupling is moderate 
below Tq, one is interested in investigating the role of 
mesonic (both thermal and quantum) fluctuations on the 
Ua{ 1) factor, using a temperature independent coupling. 

The goal of this paper is to develop a new method that 
is able to give an account of the effect of the mesonic 
fluctuations on the anomaly. This will be achieved by 
generalizing the chiral invariant expansion 0 of the ef¬ 
fective potential in the functional renormalization group 
(FRG) formalism. 

FRG and related approximations have been proven to 
be valuable tools in understanding the phase structure 
of scalar theories prl - Eoj . It has been shown in several 
studies that, even the leading order of the derivative ex¬ 
pansion applied at finite temperature gives decent results 
21-l27|. which makes the approach convenient and reli¬ 
able. The chiral invariant expansion technique, devel¬ 
oped in 0 , has the numerical advantage that instead of 
calculating the effective potential in a three-dimensional 
grid, one may derive flow equations for one-dinrensional 
coefficient functions. Furthermore, it should be noted 
that the FRG technique has not yet been applied widely 
to theories with nonvanishing Ua{ 1) anomaly for Nf = 3. 
A first analysis was presented in [28| for gauge theories, 
followed by 0] for the linear sigma model, also with the 
inclusion of quark degrees of freedom (30]. These studies 
provided important results on the effect of the anomaly 
factor; nevertheless, the flow of the anomaly coefficient 
and its temperature dependence in scalar theories have 
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not been considered in the literature before. In this pa¬ 
per we are developing a method that besides providing 
the finite temperature flow of the ‘t Hooft term also re¬ 
alizes a resummation of a wide class of U A ( 1) breaking 
operators, as will be explained later. We emphasize that 
building phenomenology upon the scheme presented here 
is beyond the scope of the paper. In this study we are 
searching for an answer for the relevance of the mesonic 
fluctuations in regard to the induced temperature depen¬ 
dence of the anomaly factor. Explicit symmetry break¬ 
ing terms representing finite quark masses are not intro¬ 
duced; nevertheless, we clarify the role of the six-point 
invariant which has not been investigated before. 

The paper is organized as follows. In Sec. II, we intro¬ 
duce the model and review the chiral invariant expansion 
technique. We go beyond the approximation presented 
in 0 , and clarify the role of the 6-point invariant of the 
theory. In Sec. Ill, we introduce the anomaly and derive 
the flow equation for the field dependent anomaly coef¬ 
ficient, which is at first analyzed analytically, and then, 
in Sec. IV the reader finds the details of its numerical 
solution. Sec. V is dedicated to conclusions. 


II. FLOW EQUATIONS AND THE CHIRAL 
INVARIANT EXPANSION 

The model to be investigated is a field theory of a 3 x 3 
matrix field <h, defined as 


$ = ( a a +iTr a ) — , 


(.a = 0...8), 


( 1 ) 


where A a are the Gell-Mann matrices, A 0 = y |1, with 
the (J a , 7 T a coefficients being scalar and pseudoscalar 
fields, respectively. $ serves as an order parameter of 
the chiral transition, and its fluctuations give account of 
the scalar and pseudoscalar meson nonets. Before we dis¬ 
cuss the details of the Ua{ 1) factor, we go through the 
anomaly free model and the approximation scheme to be 
employed. The Lagrangian is 

£ = 9' i $ t 9 /i $ _ m 2 Tr ($1$) 

- Ai[Tr($ t $)] 2 - AzTr^iW-I-), (2) 

which is clearly invariant under chiral U( 3) x 17(3) trans¬ 
formations. We choose m 2 < 0 and A2 > 0 (and also 
Ai + 3A 2 > 0), being necessary conditions that lead to 
the expected symmetry breaking pattern 17(3) x U (3) —> 
17(3), realized by <I> = uoA°/2; see details in pill ]. 

For obtaining the effective potential of the model, 
we employ the functional renormalization group method 
17,|l9(. The central object in the formalism is the scale 
dependent effective action T*,, which includes fluctuations 
with momenta q > k, obeying the flow equation 


d k T k = - Tr 


(Ti 2) + R k )~ l d k R k 


(3) 


( 2 ) 

where Tj,' is the second functional derivative of T k with 
respect to $, and R k is an appropriately chosen regula¬ 
tor function. It can be easily shown that if the UV cutoff 
is denoted by A, the scale dependent T k functional in¬ 
terpolates between the classical action [k = A] and the 
quantum effective action [k = 0]. 

In this study we use the local potential approximation, 
which is the leading order contribution of the derivative 
expansion; in other words the effective action at all scales 
is approximated as 

T fc [$] = J d 4 x (d^(x)'d^(x) - V k (x ; $)) , (4) 

where 14 is called the scale dependent effective potential. 
We use Litim’s regulator ||32| : 

R k (po,p) = (k 2 - p 2 )0(fc 2 - p 2 ), (5) 


which, at finite temperature T, leads to 


d k V k [<f>] 


67T 2 


t EE 


k 2 


m' 


( 6 ) 


There is a finite temperature sum over bosonic Matsub- 
ara frequencies uij = 2njT , and another one correspond¬ 
ing to the excitation spectrum. Here m 2 {k) denotes the 
eigenvalues of the scalar and pseudoscalar mass matrices 
= d 2 V k /d<7 l d<j J and m 2 ^(k) = d 2 V k /dn l dn : ’, 
respectively. 

Now we review the chiral invariant expansion technique 
developed in 0- The most important observation is 
that the V k local potential must reflect the 17(3) x 17(3) 
symmetry of the theory, which means that its variables 
are actually group invariants: 


V k @) = V k (h,h,h), (7) 


where 


I x = Tr (47$), I 2 = Tr [<F'4> - Tr ($ t $)/3] 2 , 

I 3 = Tr [47$ - Tr («> t 3>)/3] 3 . (8) 

Since vector symmetries cannot be broken spontaneously 
[33J, we expect $ ~ 1 and hence only I\ takes nonzero 
value in the vacuum. Motivated by this observation, one 
attempts to expand V k (Ii,l 2 ,h) around T4(/i,0,0) = 
U(I\) to get 

V k (h,I 2 ,h) « U k {h)+C k (h)h+D k {h)h + (9) 

where C k {h) = dI4/<9/ 2 | j 2 =j 3 =o and D k {h) = 

914/9/31/ 2= / 3= o. At the UV scale A, in accordance with 
©, we identify the coefficient functions with the follow¬ 
ing combination of the bare coupling constants: 

U A = mlh + (A 1A + E^)/?, C A (h) = A 2A , (10) 

with D a = 0. Higher order contributions in expansion 
© are going to be neglected, but note that, this already 
goes beyond the previous attempt Our task is to de¬ 
rive flow equations for U k (Ii ), C k (li) and D k (Ii). Note 
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m 2 /A 2 = -0.05, X, = -0.0889, X 2 = 0.4667 m 2 /A 2 = -0.05, X 1 = -0.0889, Xg = 0.4667 




FIG. 1. Demonstration of the effect of the third invariant (I 3 ), as a function of the temperature. The plots show that I 3 has 
no effect on the effective potential, and only a moderate effect on the next-to-leading (NLO) order coefficient Ck for small vo 
values. The temperature is varied between T /A = 0 — 0.6 and increases on the left (decreases on the right) from bottom to top. 


that, since the mass eigenvalues [m 2 (k)] are not chiral 
invariants, it is not straightforward to make © compati¬ 
ble with (pTJ) . A detailed description can be found in 0; 
here we just sketch the procedure shortly. 

As already stressed, in the case of the symmetry break¬ 
ing realized by $ = vq \°/2 ~ 1, J 2 and I 3 vanishes; there¬ 
fore, if we are to derive flow equations for the respective 
coefficient functions [i.e., Ck{h) and Dk{Ii)], we have to 
assume the existence of a more general condensate. We 
choose to include an infinitesimal piece proportional to 
the matrix A 8 : 


$ = u 0 A°/2 + u 8 A 8 /2. (11) 

Since 

= (12a) 

h\ vo ,v s =^(v & - 2 V 2 v 0 ) 2 , (12b) 

h\v o , V 8 = 4 ^(v 8 - 2 V 2 v 0 ) 3 , (12c) 


the flow of each coefficient function U k (h), C k {Ii), and 
D k {I\) can be easily identified; we just have to calculate 
the mass spectrum in the background m (see the ap¬ 
pendix) , substitute it into the right-hand side of © , and 
expand it around v$ = 0. Based on ©, the 0(1) term 
will give dkUk, O(uf) leads to the identification of I 2 and 
therefore dkCk, and finally O(uf) provide I 3 and d k D k - 


These read as follows: 


dkUk(h) 


dkC k (h) 


k A Ty, 
67t 2 •4—' 

U)j 

9 8 

u 2 + El u, 2 + El 

k 4 Ty, 
6 n 2 2 —/ 

Ldj 

’ 4(3Cfc + 2IiC' k ) 2 /3 

{u 2 + Elf{u 2 +ED 


u 2 j+El J’ 

(13a) 


128C7|/?/3 24 C k {C k -IiC' k ) 

+ ( W ? + E 2 f(ui 2 + El) 3 + ( W f + El) 3 
4 (3 C k C' k h + M\C' k + C k (3C k - 2 C''/ 2 )) /3 
( u * + El)(u,* + E*) 2 
MCplCk-hCl)/3 48 C 2 I 2 Cj 

K 2 + ^) 2 K 2 + ^ 0 ) 3 H 2 + El){^ + Elf 

, 6C fc - \7I\C' k 1 6 C k + 9/1 C' k + 2 IlC'l 1 

+ (^+Elf h W + Elf h 

4C'fc(6Cfc + 9 hC' k + 2I 1 2 C'(()/3l 

(^ + El)(^ + Elf J’ 1 


where a , ao and 7 r denote the scalar and pseudoscalar 
excitation spectrum, belonging to the breaking (7(3) x 
(7(3) —1 (7(3). The energies are 

E% = k 2 + M 2 = k 2 + U' k {h), (14a) 

El o =k 2 + M 2 ao =k 2 + U' k (h) + bi C k (h), (14b) 

E 2 = k 2 + M 2 a = k 2 + UKh) + 2I 1 UZ(I 1 ), (14c) 


with multiplicities 9 ® 8 ® 1, respectively. The equation 
of d k D k is too lengthy and we do not list it here ex- 
plicity; it can be found in the appendix. All Matsubara 
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m 2 /A 2 = -0.05, = -0.0889, A 2 = 0.4667 



V A 


m 2 /A 2 = -0.05, >., = -0.0889, ;. 2 = 0.4667 



FIG. 2. Comparison of next-to-leading (NLO) and next-to-next-to-leading (NNLO) coefficients of the chiral invariant expansion 
at zero temperature. Numerical values show that Dk 1 1 Ck, which means that the effect of the six-point invariant is negligible 
compared to that of the four-point invariant. 


sums appearing in the coupled flow equations can be per¬ 
formed analytically with the corresponding formulas also 
presented in the Appendix. 

An important issue of the invariant expansion ([2]) is 
the investigation of its stability, which was not carried 
out in 0- Going beyond this earlier work, i.e. the 
investigation of the flow equations of Uk and Ck, now 
we are in a position to clarify the relevance of Dk- A 
few lines on the numerics can be found in Sec. V; here 
we just shortly review the results on the stability of the 
chiral invariant expansion. 

In Fig. HI we plot a typical solution of the tempera¬ 
ture dependence of the effective potential at k = 0.2A, 
with and without the six-point invariant term included; 
in other words we investigate how the solutions change 
if we drop the equation for Dk . One observes that the 
points are on top of each other, showing that the six-point 
invariant do not play any important role. This can be un¬ 
derstood from Fig. [2 where we compare the relevance of 
Dk with Ck- Based on dimensional grounds, if one is in¬ 
terested in the effects of the NLO (Ck) and NNLO (Dk) 
coefficients on the effective potential, one actually needs 
to compare Dk • I\ with Ck ■ The figure shows that since 
at the UV scale, D^a is initiated as zero, its flow does 
not lead to a significant increase. Note that, this behav¬ 
ior is nontrivial, since there is no IR stable fixed point in 
the theory that would lead to scaling and (trivially) the 
suppression of nonrenormalizable couplings in the IR. 

We expect that since the six-point invariant does not 
play a crucial role, none of the higher order Taylor co¬ 
efficients have any effect on the solution either. With 
these findings the chiral invariant expansion has proven 
to be stable, and one can safely truncate the series at 
next-to-leading order. 


III. INCLUSION OF THE U A ( 1) FACTOR 

Implementation of the U A (V) anomaly is done via ‘t 
Hooft’s determinant term. We add the following term 
into Lagrangian: 

£-UaO = a ' (det + det 4>), (15) 

which explicitly breaks the Ua( 1) subgroup of the U (3) x 
U( 3) chiral symmetry. The anomaly changes the spec¬ 
trum to8©l®7®l®l [see also (IA2D in the appendix], 
and the flow equation ([H]) generates all Ua( 1) breaking 
operators, such as, e.g., Tr (<J>t<f>) (det <I>t + det 4>) into 
the effective potential. 

The effect of the anomaly on the effective potential can 
be formulated as 

Vk = I4|a=o(-fi) h, h) 

+ (h, h, / 3 )(det & + det $)\ (16) 

i 

where 14-| a =o is the anomaly free effective potential. 
Since the anomaly can only be carried by the /det = 
det$t +det$ operator, the coefficients are U( 3) x 
U( 3) invariants, and thus one can apply the chiral invari¬ 
ant expansion on them as well: 

A®(I 1 ,I 2 ,I 3 ) = A®(I 1 , 0,0)+ £ n A kU^K 3 - 

{a} f—2,3 

(17) 

In this paper we only present a leading order analysis: 

Vk w T4| a =o + A A) (h,I 2 ,h)( det & + det $), (18) 
and furthermore, we neglect the / 2 and I 3 dependence of 
Ak, 

A A) (h,h,h) * a£\i u 0 , 0 ) = Ak(h). (19) 
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Note that a field dependent Ua( 1) coefficient already 
goes beyond a simple perturbative renormalization group 
analysis of the anomaly. Taylor series of A k (I 1 ) around 
I\ = 0 shows that an infinite resummation of operators 
[TV (<fd<£>)]"(det $1 + det <f>) realizes in the system. 

As already announced in the introduction, in the lit¬ 
erature the anomaly coefficient (a) is somewhat contro¬ 
versially linked to the topological susceptibility. Here we 
treat it as a temperature independent parameter and in¬ 
vestigate how thermal fluctuations of mesons affect the 
strength of A k {lY). 

In order to obtain an approximate flow equation com¬ 
patible with GHD-dUD, we have to expand J6|) around the 
zero anomaly configuration, i.e., A k (I\) = 0. The reader 
is referred to the appendix for formulas of the mass ma¬ 
trices and the necessary derivatives of invariants. The 
zeroth order of the expansion gives the already obtained 
flow equations for U k , C k , and D k , while the next-to- 
leading order terms combine into a term proportional to 
/det- The coefficient of /det becomes the scale derivative 
of Afc(/i). We arrive at 


d k Mh) 


T 


\ ' 


9 M 


9 A k 


K + 


67T 2 2 —-' 

8 A' k 

W + E* Y ' MuME *) 2 


+ 


KY 

12 A k 


H^ + eiy 


w 3 1 a o> 

3 A k 


+ 


7 A' 


+ 


2 hK 


W + ElY-h {u£+E%Y W+El) 


( 20 ) 


This equation can now be solved numerically given that 
the functions U k , C k and D k are known. 

Before we present the full numerical results, let us solve 
(Blill using the assumption of V k taking the form of the 
classical (bare) potential, with /-dependent, field inde¬ 
pendent couplings. For constant A k , at high enough tem¬ 
peratures, around I\ ~ 0 (f20li simplifies: 


8 k log A k = 


2k 4 T 8 C k - 3 U” 
3^ 2 (/ 2 + t /') 3 ' 


( 21 ) 


The effective potential in our approximation, V k = 
U k (h) + C k (h) ■ h + D k (h) ■ h + A k {I\) ■ /det can be 
parametrized as [see also USD] 

U k (h) = mlh + (Ai fc + ^f)ll 

C k (h) = X'2k, D k (Ii) = 0, 

A k {I\) = a k , (22) 


from which at high temperature 


®k — 7 


L dK 


(^ + U ' K )3 ( A 2*-Ai* 


(23) 


This result shows that within the field independent cou¬ 
pling approximation, and with the assumption that the 
flow is monotonic, the condition of decreasing anomaly 


strength as a function of the temperature can be formu¬ 
lated already at bare level: 

A 2 A < Aia- (24) 


IV. NUMERICAL RESULTS 

Equation (1231) might not survive exactly in the numer¬ 
ics, but it turns out to be a qualitatively correct approxi¬ 
mation. The applied method for solving the coupled flow 
equations (ll3aD - (ll3bl) . (IBID , and (EOl) has been the same 
as in M- At all scales (fc) the U k (h), C k (I \), D k (Ii ), 
and A k (Ii) functions are stored on a grid, typically in 
the region [0,2], with a step size of ~ 1CU 3 , and all di¬ 
mensionful units are measured in terms of the UV cutoff 
A. The flow equations are integrated with the Runge- 
Kutta method with a typical step size of Ak ~ 10 -5 . 
We also emphasize that, due to numerical stability the 
flows have to be stopped [l(| around the scale where the 
potential is g radually becoming convex [ij|, |34j, [35]. As 
shown in |10j], critical temperatures and the correspond¬ 
ing discontinuities of the order parameter (in case of first 
order transitions) can be obtained by extrapolation to 
k = 0 from the k > 0 results. Note that, in the present 
study however, we decided to stop the flows uniformly at 
k = 0.2A, since we are only interested in the tendency 
of the evolution of the temperature and field dependent 
anomaly, which can perfectly be seen around k ~ 0.2A. 
Extrapolated results to k = 0 with physical parametriza- 
tion of the model will be reported elsewhere. 

In Fig. H we show the £Ai(l) coefficient function at 
two distinct points of the parameter space that lead to 
either anomaly strengthening or weaking. We observe 
that already at zero temperature, quantum fluctuations 
develop a structure for A k (Ii) showing that a field in¬ 
dependent approximation is not appropriate. It can also 
be seen that the full numerical solution is in accordance 
with (J2U); at I\ ~ 0, for large A 2 the anomaly is getting 
weaker as a function of the temperature, while if Ai dom¬ 
inates, it is strengthening. It is interesting to see that for 
larger I\ values the former statement is not necessarily 
true anymore: even though at small values of the field the 
anomaly weakens, if I\ is larger, it might also strengthen 
at the same time. In case of large A 2 , we observe that 
there might be a maximum value for the ZA(1) factor as 
a function of the temperature, followed by a moderate 
decrease. 

The vacuum state is of course always determined by 
the effective potential; therefore, in Fig. [3] we also in¬ 
dicate the field values that minimize the complete effec¬ 
tive potential (i.e., V k ). This shows that a field indepen¬ 
dent anomaly coefficient [belonging to the leading order 
of the Taylor expansion of A k (Ii) around I\ ~ 0] is in¬ 
deed a crude approximation; at some temperatures the 
actual strength of the anomaly in the minimum may dif¬ 
fer up to even ~50% from the value of it in the origin. 
Based on this observation, several possibilities can arise 
as the temperature is raised. The anomaly can either 
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m 2 /A 2 = -0.05, X-, = 0.1667, ^ = 0.1 m 2 /A 2 = -0.05, A-, = -0.0889, ^ = 0.4667 




V 0 /A V 0 /A 


FIG. 3. Temperature dependence of the anomaly coefficient at k — 0.2A. The plots show that the rough analytical estimate 
is qualitatively correct, depending on the sign of (A 2 A — Aia) the anomaly can either weaken or strengthen due to mesonic 
fluctuations. Horizontal dotted lines represent the anomaly coefficient at the UV scale, and the circle-shaped spots show the 
field values that minimize the effective potential, together with the corresponding anomaly strengths. 


weaken or strengthen (more or less) nronotonically, or it 
can increase and discontinuously drop at the transition 
point. The latter can be snroothened by explicit symme¬ 
try breaking terms representing finite quark masses, or 
by the initial strength of the anomaly itself. The map¬ 
ping of this issue will be reported elsewhere. 

We note that recent tree-level parametrizations of the 
model (36j, or extensions of the model [37, HU, show that 
A 2 > Ai. It is important to stress that one should not 
draw any final conclusions on the behavior of the anomaly 
based on these findings, since the value of the bare cou¬ 
pling constants strongly depend on the cutoff scale and 
regularization, and also on the employed approximations. 
There are other recent parametrizations of the model, us¬ 
ing the FRG formalism, where Ai ~ A 2 [dll ■ 

Finally let us point out that since the I\ invariant spec¬ 
ifies an 0(18) combination of the fields [7i ~ (<r“) 2 + 
(7r“) 2 ], roughly speaking Ai characterizes the orthogonal 
nature of the theory, while A 2 shows its deviation from 
it. Based on this argument, one can reformulate (1M|) as 
a requirement of the theory to be closer to the 0(18) 
model, if one is to seek decreasing anomaly factor as a 
function of the temperature. 


V. CONCLUSIONS 

In this study we have analyzed the three flavor lin¬ 
ear sigma model with the local potential approximation 
of the functional renormalization group in 3 + 1 dimen¬ 


sions. The chiral invariant expansion introduced in M 
has been developed to the case of nonzero Ua{ 1) anomaly, 
and its stability has also been checked explicitly. At 
first, by deriving the flow equation of the next-to-next- 
to-leading coefficient of the expansion we have found that 
the chiral invariant expansion does appear to be stable; 
there is no need to go beyond next-to-leading order. Fur¬ 
thermore, we derived the flow equation of the field depen¬ 
dent Ua( 1) anomaly factor, and found qualitatively that 
the first quartic coupling has to dominate the second for 
decreasing anomaly as a function of the temperature. 

It has turned out that mesonic fluctuations are capa¬ 
ble of describing anomaly weakening, even without an 
implicit temperature dependence of the ‘t Hooft deter¬ 
minant coupling. Adding explicit symmetry breaking 
terms, and parametrizing the model with low energy 
mesonic spectrum would allow us to investigate temper¬ 
ature dependence of the anomaly and spectrum itself in 
the physical point. This represents a future study to be 
reported in the near future. 
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Appendix A. MASS MATRICES AND GROUP INVARIANTS 


In this appendix we give all the necessary formulas to calculate the scalar and pseudoscalar mass matrices m 2 a l3 (k) = 
and m 2 i j(k) = They are needed to obtain the flow equations for £4, C k and Dk , and ultimately for 

Ak- 


dl\ dl\ 
d & 1 doi 
dli dh 
da 1 d erf 
dh dh 


• h+mo ■ 6 ) + £qh(Wi)+mo ■ h+ mo ■ h) 


dh dl 2 
dai da 1 


1)C1(A) 


(dhdh_ 

V da * 9crf 




( 

I^IS-( c/i ' <jri,+c ' i ' <Jl) ' h + mo ■ * 


da l dai 

dh dh 
dai da 1 
d 2 h 


)mo, 




+ 


d^diri 


da‘dai da l dai 

(y'k(h) + C' k (h) ■ h + D' k (h) ■ 7 3 ) 


D k (h), (Ala) 


+ 


dh dh_ dh dh 
di r* dni diri dn 1 


)rno +(grU+HH) AU.) + Stimo + Stjmo- m) 


d 2 h 


If we include the ‘t Hooft determinant, and the corresponding A k (h) ■ (det <I> + det<f>f) = A k {h) • Idet term in the 
scale dependent effective potential, the mass matrices get corrected by 


Am cr,ij g a i Q a j A k (h ) ■ Idet + 


d 2 h 






A'k (II ) • Idet + ( 

di: l d^ Ak ^ ■ /det + (dir* dir* 


da l dai 

d 2 h 


dh dlfe t dh dhe t 

da 1 da j da ] da 

dh dl det , ( 9 /l dldet 


det 


-)a , / s ( Ji ) + 

dni dn 1 ) dir 1 dir* 


d 2 I t 
da l dai 
d 2 Idet 


A k (h),( A2a) 
A fe (J 1 ).(A2b) 


Invariants /i and 4, and their first derivatives in the two-component background are the following: 


vl + vi 


L 1 bo,«8 


hlvotVs — 9^(^ 8 2y/2vo) , 


dh 

da a 

dh 

da a 


= v 0 S aU + v 8 S aS , 


Vo,Vs 


( 2 v 0 v £ 

V 3 2 ,V 2 




dh 

dn a 

dh 

(?7r a 


= 0, 


= 0. 


The second derivatives are 

d 2 h 


da a da b 


= S' 


ab 


d 2 h 


dTr a dir b 


= 6 ' 


ab 




ff w 8» 



if 

a = 

= b 

= 0 




V2 + 3 

V 0 V 8 , 

if 

a = 

= 0 , 

6 = 8 

or a 

d 2 I 2 

— > 

2,,2 4- 
1 3^0 + 

«8 

2 

- V2v 0 vs, 

if 

a = 

-b- 

= 8 


da a a b 

— S 

^0,^8 

2„,2 , 
3^0 + 

£ 

6 

+ VZvovs, 

if 

a = 

-b- 

= 1,2, 

3 



2 2 , 
3^0 + 

6 

- j^v 0 v 8 , 

if 

a = 

= b 

= 4,5, 

6,7 



,0, 


else 






'0, 



if 

a = 

= b 

= 0 




vl 

— 3 V 2 

+ j 

5 V 0 VS, 

if 

a = 

= 0 . 

6 = 8 

or a 

d 2 h 

= < 

I ^ - ^0^8, 

if 

a - 

= b 

= 8 


dTT a n b 

Vq,V 8 

-£ + 
6 ^ 

^V 0 v 8 , 

if 

a = 

= b 

= 1,2, 

3 



1^8 - 

3j2 V ° V »’ 

if 

a - 

= b 

= 4,5, 

6,7 



, 0 , 



else. 





(A3a) 

(A3b) 

(A3c) 


(A4a) 


(A4b) 


(A4c) 



Furthermore, invariants I 8 and Idet, and their first derivatives are 


^3 |u 0 ,^8 ~ 288 2 "\/ 2 no) , 


Id*!™ = A - - A(A5a) 


dh 

da a 

dldet 


Vo, v 8 


dcr a 


.)> _ | ( „ 8 _ 2V2vof(V2va 
2V 0 ~ V {§aO _ MV2vo + V 8 ) 5 g S 


§ 
dldet 


2 y /6 


Vo,Vs 

Finally, the second derivatives are 

,,3 


2 V 3 


d 2 h 


diT a a b 


= 


d 2 h 


dTT a TT b 




^/g( v 8 - 2 t/2v 0 ), 

-§(24s/2v$ - 32v 0 V8 + 5 V2v§), 

H(—16-\/2uo + 48^0^8 - 20v / 2vqV8 + 5uf), 
H(16v / 2^o + 16 vfiv 8 + 4\Z2v 0 v§ - 5v§), 
|§(-8v / 2co + 16v$v 8 - 5y/2v 0 v§ + uf), 

vO, 

0 , 

_ 4t/2v 0 v 8 + i’l), 
i| ( 8 v§ - 4y/2vov 8 +vD, 

-3|( 8u o - 4\/2vov 8 + vD, 

2§(4wo “ 11\/2uov 8 + 5vg), 

0 , 


d 2 L 


det 


da a a b 


\J 3 V °’ 


_V8_ 


V6’ 


Vo 

118 

V6 

x/3’ 

Vo 1 

118 

V6 + 

%/S’ 

Vo 

118 

V3 

2^3 


v«, 


d 2 I t 


det 


dTT a TT b 


y 3 U 0; 

VS 

Ve’ 

Vo I _«8. 

= V6 V3’ 

V0_ _ V8_ 
VE y/3’ 

VO 1 1>8 

^ 2 V 3 ’ 

0. 


3 ^ 2-\/6 6 a /3 

= 0, (A5b) 


dit a 


= 0 . 


if a - 

if a ■ 
if a - 
if a ■ 
if a 
else 

if a - 
if a - 
if a - 
if a - 
if a - 
else 

if a = 
if a = 
if a = 
if a = 
if a = 
else 

if a = 
if a = 
if a = 
if a = 
if a = 
else 


= 6 = 0 

= 0, 6 = 8 or a = 8, 6 = 0 
= 6 = 8 
= 6= 1,2,3 
= 6 = 4,5,6,7 

= 6 = 0 

= 0, 6 = 8 or a = 8, 6 = 0 
= 6 = 8 
= 6= 1,2,3 
= 6 = 4,5,6,7 

6 = 0 

: 0, 6 = 8 or a = 8, 6=0 
: 6 = 8 
= 6= 1,2,3 
6 = 4,5,6, 7 

:6 = 0 

0, 6 = 8 or a = 8, 6 = 0 
6 = 8 
6= 1,2,3 
6 = 4,5,6,7 


(A5c) 


(A6a) 


(A6b) 


(A6c) 


(A6d) 


Appendix B. NNLO FLOW EQUATION 

Here we present the flow of the coefficient function Dk(Ii); for explanation see ©, J9J, and (fl4l) . 


a ‘ D * = 


Cl 


C2 


C3 


C4 


K 2 + ElY 


(u 2 + ED 2 (u 2 +El) (rf + E?) 


C 5 


C6 


C7 


^ 1 i^+Ein^+ED 

C 8 


0 '^ + El 0 ) 3 K 2 + ^ 0 ) 4 i^+E^n^ + ED (u ; 2 + ED 2 

eg , cio , Cn 


+ 


C 12 


(^ + M2)3( W 2 + M 2 )2 ^ + El) 3 ( W ?+£3) 3 (u,? + £S ) 3 (a ,2 + £2)4 


(Bl) 
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where the Cj coefficients are the following: 

Cl =(12C 2 D k h - 24 C k D k I 2 C' k - 48 D k I*C? + 20 C 2 I 2 D' k + 48 l\C k C' k D' k - 8 C 2 /-?!?" - 243C 2 t/" 
-90 G k D k hU'l - 324 hC k C' k U£ - l 8 QI 2 D k C' k U'y - 108/ 2 C' 2 [/" + m\C k D' k U'; - 72 llC' k D' k U'l 


+24lfC k D k U k - 135 I k D k U' k * - 117/ 2 D'[/f - lSJ 3 4 !)^)// 2 ^ - 3[/£) 2 , (B2a) 

c 2 =-18Cf//i - 24C k C k - 8/iCf, (B2b) 

c 3 ={WC 2 k D k h + 21QC k IiC' k + 9 QC k D k I 2 C' k + 192 D k lfC[ 2 - 416 C 2 k I 2 D' k - 192 C k l\C' k D' k + 972Cf C/j" 
+mC k D k hU k + 648IiC k C k U k + 720D k lfC k U k + 432/ 2 Gft/" + miC k I 2 D’ k U{ + 288lfC' k D' k U £' 
+788D k I 1 U k 2 + 486/iCfcl7fc 2 - 288/ 1 2 D , fc 17fc 2 )/4/ 1 2 (2C' fc - 3U k ) 2 , (B2c) 

c 4 =-18C^//i - 24C k C k - 8hC' k 2 (B2d) 

c 5 =(18C 3 - 4 C k D 2 I 2 + 180 hC 2 C' k + 96 1 2 D k C k C' k + 48I 2 C k C' 2 + 32 D k C' k l\ + 135C 2 U" 

+108 C k D k hU” + 6D k I 2 U k - MhCkC'kU'^/h(C k - 3U"/2), (B2e) 

c 6 =2(27 Cl + MC 2 k D k h + 36 C k D 2 I 2 + 8 D 3 k lf), (B2f) 

c 7 =72 hC 2 C' k , (B2g) 

c 8 =-27£> fc /4Ji - 27C' k /2I\ - 9 D' k , (B2h) 

c 9 =—32C k Ii + 32I 2 C k C k , (B2i) 

c w =18C k (D k + C k /h), (B2j) 

on =—64C k I 2 /3, (B2k) 

ci2 =-6 Cl (B21) 


Appendix C. MATSUBARA SUMS 


The Matsubara sums appearing in the flow equations have the following form: 

1 


S(i,j) = 


(u,l + Efy(u,l+Ei)r 


First, we calculate 5(1,0) and 5(1,1): 


5(1,0) 

5(1,1) 


coth(£i/2T) 

2E\ ’ 

1 Ex coth(£’ 2 /2T) - E 2 coth(£i/2T) 
2ExE 2 E 2 - El 


The rest can be obtained by differentiation: 

S(n > 1, 0) 
S(n > 1, to > 1) 


(-*)”“ 1 0 " _ 1 5 ( 1 , O ) 

(n-1)! d(E 2 J"- 1 ’ 

(-*)”- 1 (- i )” 1 " 1 d "- 1 ^“^( l , 1 ) 
(n — 1)! (to — 1)! d(Ef) n ^ 1 d(E 2 ) m ~ 1 


(Cl) 

(C2) 

(C3) 

(C4) 

(C5) 
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